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ABSTRACT 

Reflected solar radiative fluxes emerging for the Earth’s top of the atmosphere are inferred 
from satellite broadband radiance measurements by applying bidirectional reflectance functions 
(BDRFs) to account for the anisotropy of the radiation field. BDRF’s are dependent upon the 
viewing geometry (i.e. solar zenith angle, view zenith angle, and relative azimuth angle), the amount 
and type of cloud cover, the condition of the intervening atmosphere, and the reflectance 
characteristics of the underlying surface. A set of operational Earth Radiation Budget Experiment 
(ERBE) BDRFs is available which was developed from the Nimbus 7 ERB (Earth Radiation Budget) 
scanner data for a three-angle grid system. An improved set of bidirectional reflectance is required 
for mission planning and data analysis of future earth radiation budget instruments, such as the 
Clouds and Earth’s Radiant Energy System (CERES), and for the enhancement of existing radiation 
budget data products. 

This study presents an analytic expression for BDRFs formulated by applying a fit to the 
ERBE operational model tabulations. A set of model coefficients applicable to any viewing condition 
is computed for an overcast and a clear sky scene over four geographical surface types: ocean, land, 
snow, and desert, and partly cloudy scenes over ocean and land. The models are smooth in terms of 
the directional angles and adhere to the principle of reciprocity, i.e., they are invariant with respect to 
the interchange of the incoming and outgoing directional angles. The analytic BDRFs and the 
radiance standard deviations are compared with the operational ERBE models and validated with 
ERBE data. The clear ocean model is validated with Dlhopolsky’s clear ocean model. Dlhopolsky 
developed a BDRF of higher angular resolution for clear sky ocean from ERBE radiances. 
Additionally, the effectiveness of the models accounting for anisotropy for various viewing 
directions is tested with the ERBE alongtract data. An area viewed from nadir and from the side give 
two different radiance measurements but should yield the same flux when converted by the BDRF. 
The analytic BDRFs are in very good qualitative agreement with the ERBE models. The overcast 
scenes exhibit constant retrieved albedo over viewing zenith angles for solar zenith angles less than 
60 degrees. The clear ocean model does not produce constant retrieved albedo over viewing zenith 
angles but gives an improvement over the ERBE operational clear sky ocean BDRF. 
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Chapter 1 
INTRODUCTION 
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In order to enhance our understanding of the radiative energy interaction 
between the Earth and space, the components of the Earth's radiation budget need 
to be examined. The availability of remotely-sensed radiance measurements from 
earth-orbiting and geostationary satellites makes this investigation feasible. The 
components of the radiation budget are the incident solar flux, the Earth s emit- 
ted (longwave) radiation, and the Earth's reflected (shortwave) solar radiation at 
the top of the atmosphere (TO A), which is considered to be the surface of reference 
(Fig. 1-1). The TO A longwave (LW) and shortwave (SW) fluxes are not directly 
measurable quantities, but rather need to be derived from the observed radiances. 
The radiation emerging from the TOA has an anisotropic distribution and is influ- 
enced by the reflectance characteristics of the underlying surface, the illumination 
and viewing conditions, the optical properties of the intervening atmosphere, and 
the amount of cloud coverage within the target area [1,2]. The radiance to TOA 
flux conversion requires knowledge of the angular characteristic of the outgoing 
radiation field described by angular distribution models (ADMs) or bidirectional 
reflectance functions (BDRFs), the latter being the term used in this study. These 
BDRFs account for the dependence of spacecraft-measured radiances on the view- 
ing geometry. Uncertainties in these models lead to errors in the derived net fluxes 
at TOA, thus describing the Earth radiation budget inaccurately. 

Extensive studies of the anisotropic reflectance characteristics for various 
surface types have been performed using instrumentation on balloons, high-alti- 
tude aircraft and on satellites. Earlier assumptions of an isotropic solar radiation, 
corresponding to a Lambertian Earth surface, resulted in significant errors in the 
computation of the reflected fluxes [3,4]. Larsen and Barkstrom [5] determined, for 
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Figure 1-1. Radiative interaction in the Earth-atmosphere system. 
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instance, that for Arctic latitudes, assumptions of an isotropic surface causes the 
global albedo to be overestimated. Coulson and Reynolds [6] used laboratory- 
gathered bidirectional reflectance data to describe the anisotropic properties for 
various land surfaces, discriminated by soil and vegetation type. Brennan and 
Bandeen [7] used a medium resolution infrared radiometer (MRIR) aboard the Na- 
tional Aeronautics and Space Administration's (NASA) Convair-990 high-altitude 
research aircraft to investigate the anisotropic patterns of clouds, water and land 
surfaces. Salomonson and Marlatt [8] studied the anisotropic characteristics of 
three highly reflective surfaces, namely stratus clouds, snow and white gypsum 
sand, from data collected using the NIMBUS F3 MRIR mounted on a Piper Twin 
Comanche aircraft with a flight altitude ceiling of 9000 meters. 

Because satellites provide more comprehensive spatial and temporal coverage 
over the entire earth, satellite measurements are ideal for developing anisotropic 
models. Moreover, since measurements are made above the atmosphere, more re- 
alistic determination of the incoming and outgoing radiation is possible. Addi- 
tionally, surfaces can be viewed from a larger range of illumination and viewing 
directions. Ruff et al. [9] studied the angular properties of clouds from radiometric 
measurements aboard the TIROS IV spacecraft. Although relatively high resolu- 
tion radiometers were deployed on the TIROS satellite series and on NIMBUS 2, 
the need to obtain data that complied with spatial resolution requirements over an 
extended period of time was not met until the deployment of the NIMBUS 3 space- 
craft. This mission provided the opportunity to study the annual global budget 
over a continuous period for the first time [2]. Bidirectional reflectance models for 
three surface types (ocean, cloud-land, snow) were applied for three solar zenith 
angle ranges. 

The first scanning instruments designed specifically to measure the radia- 
tion budget flew aboard the NIMBUS 6 and 7 satellites [10,11]. The NIMBUS 7 
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satellite was launched in October 1978 into a noon (ascending node) sunsynchro- 
nous orbit at an altitude of 950 km. The NIMBUS 7 Earth Radiation Budget (ERB) 
scanning instrument was equipped with a biaxially scanning set of eight optical 
telescopes, four of which monitored broadband SW (0.2 - 4.0 pm) radiation and the 
remaining four monitoring broadband LW (4.0 - 50 pm) radiation. The spatial res- 
olution varied from 90 km x 90 km at nadir to 250 km x 250 km at a maximum scan 
angle of 72* [11]. The operational scan modes of the radiometers permitted radi- 
ance observations in the principal plane of the sun where the angular variations of 
the reflected sunlight are most prominent. This instrument was designed to yield 
data necessary for the development of a comprehensive set of BDRFs. However, 
orbit configuration and scanning constraints have inhibited a more complete an- 
gular coverage. Thescanninginstrumentson the NIMBUS 6 and 7 satellites have 
provided extensive datasets from which anisotropic models have been developed. 
Taylor and Stowe [12] developed bidirectional models for eight uniform surface 
types from the NIMBUS 7 measurements. Staylor [13] and Staylor and Suttles [14] 
developed BDRFs for clouds and deserts, respectively, using NIMBUS 7 scanner 
measurements. The present models used to process data from the Earth Radiation 
Budget Experiment (ERBE) [15] were derived from the NIMBUS 7 Earth Radiation 
Budget (ERB) and the Geostationary Operational Environmental Satellite (GOES) 
data sets [16]. 

Some uncertainties in early earth radiation budget measurements are attrib- 
uted to poor temporal coverage by a single, sunsynchronous spacecraft. The issue 
of poor diurnal sampling was addressed with the inception of ERBE, a multiple 
satellite mission. The orbital configuration of the multi-satellite system allowed 
for adequate temporal and spatial coverage. ERBE consists of three spacecrafts, 
the Earth Radiation Budget Satellite (ERBS), and the NOAA-9 and NOAA-10 
spacecrafts. ERBS is a ERBE-dedicated spacecraft that was launched into orbit in 
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October 1984 in a 57' inclination orbit with an altitude of 600 km, restricting latitu- 
dinal coverage between 57*N and 57*S. The NOAA-9 and NOAA-10 are high alti- 
tude sunsynchronous satellites deployed to altitudes of 812 km and 830 km, 
respectively. The ERBS has a westward precessional period of 72 days around the 
earth while the NOAA-9 and NOAA-10 have equator crossing times of 1430 (as- 
cending node) and 0730 (descending node), respectively. The ERBE scanning 
package consists of a shortwave (0.2 - 5.0 pm), a longwave (5 - 50 pm), and a total 
(0.2 - 50 pm) radiometer. The radiometers, sweeping from horizon to horizon, are 
mounted on a single scan head and boresighted [17] so as to detect radiation from 
the same field of view. The normal operating mode is in the cross-track scan direc- 
tion; although, for limited periods in 1985 (January and August), the scanner was 
rotated 90* in azimuth to operate in the alongtrack scan mode. Operation in this 
mode is ideal for determining LW limb-darkening functions [18,19] and for vali- 
dating shortwave anisotropic models because it allows for a single site on the sat- 
ellite ground track to be viewed from different viewing zenith angles during a 
single orbital pass. 

Dlhopolsky [20] utilized ERBE shortwave radiance measurements to gener- 
ate an improved set of angular directional models with increased angular resolu- 
tion for clear sky over ocean surface. 

The accuracy of future radiation budget data products can be enhanced 
with an improved set of BDRFs. The improved set of BDRFs must be continuous 
and smooth from one angular bin to another. It must also satisfy reciprocity (i.e. 
interchanging the incident and reflected directions must yield the same flux con- 
tribution). Additionally, radiances that are measured from different viewing an- 
gles over a single site must be converted to the same flux if the bidirectional 
reflectance function is modeled correctly (i.e. there is no albedo growth from nadir 
to limb). The present ERBE operational models are discontinuous from one bin to 
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another and do not satisfy reciprocity. Moreover, there is about a 10% albedo 
growth from nadir to limb. 

The objective of this study is to develop an analytic BDRF that shows the 
dependence of the solar reflected radiation on surface, cloud, viewing geometry, 
and atmospheric conditions. This study presents a simple analytic formulation of 
the BDRFs to model the anisotropy of nine of the twelve ERBE scene types. The 
nine basic scene types are classified according to geographical surface type (ocean, 
land, snow, desert) under varying degrees of cloud cover (clear, partly cloudy, 
mostly cloudy, overcast). The remaining three scenes are mixed scenes that are as- 
sumed to be made up of 50% ocean and 50% land (e.g. coast). The model coeffi- 
cients are derived by applying an analytic fit to the NIMBUS 7 BDRFs [16] which 
were used to process ERBE data. For each scene type, a single set of model param- 
eters is required for application to any combination of viewing geometries. This 
analytic formulation satisfies the principle of reciprocity and avoids the disconti- 
nuity from one discrete angular bin to another as observed from the ERBE opera- 
tional BDRFs. The analytic BDRFs are validated by comparison of the resulting 
radiances and TOA fluxes with ERBE observations. Model results and validation 
thereof are presented for clear, partly cloudy, mostly cloudy over ocean and over- 
cast scenes. Results of this study will be used for mission-planning and data inter- 
pretation of next-generation earth radiation budget programs such as the Clouds 
and Earth's Radiant Energy System (CERES) mission [21]. 

The physical problem and theoretical formulation are presented in Chap. 2. 
The analytic form of the BDRF is developed in Chap. 3. Chapter 4 discusses the 
validation techniques and results. Finally, conclusions that were drawn from this 
study are presented in Chap. 5. 



Chapter 2 

PHYSICAL PROBLEM AND THEORETICAL FORMULATION 
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2-1. Surface Anisotropy 

When solar radiation impinges on the Earth-atmosphere system, it is reflected 
in various directions. The reflected shortwave radiances are dependent upon the 
direction from which a surface is being viewed as well as on the surface angular 
reflectance characteristics. Most natural surfaces exhibit varying degrees of reflec- 
tance in different directions. A specular surface, such as the ocean surface, is mir- 
ror-like in its nature of reflectance and highly direction-dependent, while a diffuse 
surface will reflect uniformly in all direction. This directional dependence of the 
radiation field is defined as anisotropy. All Earth surfaces exhibit some degree of 
reflectance anisotropy which is described by bidirectional reflectance functions 
(BDRFs). Aside from their dependence on the underlying surface, the BDRFs are 
also influenced by the amount of cloud cover and the state of the intervening 
atmosphere. 


2-2. Bidirectional Reflectance Model 

The reflected radiance L leaving the top of the Earth's atmosphere varies with 
direction. The geometric relationship used in this study is depicted in Fig. 2-1 
[16]. The angle between ray from the sun and the normal to the target area is the 
solar zenith angle £ while the angle between zenith ray and the normal to the tar- 
get area is the viewing zenith angle 0 . The relative azimuth angle (j) is the angular 
distance of the satellite from the principal plane, i.e. the plane containing the sun, 
the Earth's center and the point of observation. The azimuth angle for an exiting 
ray is measured from the principal plane on the side away from the sun. Thus, 
reflection in the forward direction corresponds to <|>=0* while backward reflection 
corresponds to <t> =180*. The reflected shortwave flux M is obtained by integrat- 
ing the radiances L over all the outgoing directions such that 
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Figure 2-1 . Geometrical relationship of the sun, satellite and target area 
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In n/2 

M= J J L(0, 0, 0 cose sin ed9d<J> (2- 1 ) 

o o 

The flux M has units of Wm' 2 while L has units of Wm' 2 sr' 1 . For an isotropic 
surface (i.e. the reflected radiation is the same in all directions), M = nL. The bidi- 
rectional reflectance function (BDRF) R , which characterizes the anisotropy of 
the reflected radiance, is the ratio of Lambertian flux obtained from the observed 
radiance in a given direction to the actual flux computed by Eq. (2.1) and is given 
by 

/?(M,0 = nL(M,0/M (2.2) 

Values of BDRF equal to unity imply a radiance measurement will provide the 
correct radiant flux assuming isotropy. Any departure from unity corresponds to 
the fractional error that is incurred if the measured radiance is used to estimate 
the radiant flux for an isotropic assumption (i.e. BDRF <1.0 (>1.0) overestimates 
(underestimates)) the flux by an amount equal to the portion lesser (greater) than 
unity). 

The normalization condition for R is derived by substituting Eq. (2.2) into Eq. 
(2.1) giving 


2n n/2 

J J /fcosSsinSdQd^ = n 
o o 

The albedo at TOA is defined as 


(2.3) 


M = Sa (0 cos C ( 2 - 4 ) 

where 5 is the solar flux, a is the albedo at TOA and C is the solar zenith angle. 

The reflected radiation adheres to the principle of reciprocity [22]. This prin- 
ciple states that for a given observation point, the positions of the spacecraft and 
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sun may be interchanged and still yield the same flux contribution. From this 
principle, it follows that 

/? (0, 4>, 0^(0 = 9)0(0) (2.5) 

The quantity r (0, <J>, £) = R (0, <J>, 0 a ( 0 is called the bidirectional reflectance. 

The BDRFs are dependent upon the underlying geographical type (ocean, 
land, snow, desert, coast), the cloud cover (i.e. clear, partly cloudy, mostly cloudy 
and overcast), and the viewing geometry (sun, spacecraft, and observation point 
geometry). Perturbations in the atmospheric conditions also strongly influence 
anisotropy (e.g. nature of scattering particles, turbidity). Figure 2-2 is a schematic 
illustration of BDRFs for a Lambertian surface and two variations which are non- 
Lambertian. Fluxes retrieved by using each of these will be different. It illus- 
trates that the use of the wrong BDRF will result in errors in retrieved fluxes. The 
directional models that are presently being used to process data from earth radia- 
tion missions, such as the ERBE, are not closely approximated by a Lambertian 
model. The dotted line is a more realistic representation of the anisotropy of the 
radiation field for most scene types in which the atmospheric perturbations that 
causes these models to deviate from the mean values are considered. This dia- 
gram shows that the set of operational models used to process earth radiation 
budget data is a more realistic representation of anisotropy than a Lambertian 
surface assumption. 

2-3. Results from Previous Anisotropy Studies 

Studies to determine the anisotropic patterns of the outgoing radiation field 
from various surface types have been performed using balloon, aircraft and satel- 
lite data. A few studies have been selected and their results presented in the fol- 
lowing sections. The first section briefly discusses results from aircraft, 
laboratory, and satellite (TIROS IV and NIMBUS 3) experiments. The second sec- 
tion describes the BDRFs generated from the NIMBUS 7 data. Finally, the third 
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Figure 2-2. Schematic of BDRFs for Lambertian and non-Lambertian surfaces. 
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section introduces the Dlhopolsky BDRF for clear ocean that were developed 
from ERBE radiances. 


2-3.1. BDRF from Laboratory, Aircraft and TIROS IV and 
NIMBUS 3 Data 

Coulson et al. [23] conducted laboratory experiments on the anisotropic 
characteristics of natural sands and soils, particularly red clay soil and white 
quartz sand. These surfaces were found to exhibit a more pronounced reflectance 
in the backward scatter direction than in the forward scatter region. No increase 
in anisotropy in the specular direction was observed. Furthermore, an increase in 
reflectance occurred as the solar zenith angle increased. 

Meanwhile, using a NIMBUS F3 medium resolution radiometer aboard a 
Piper Twin Comanche aircraft, Salomonson and Marlatt [8] determined the aniso- 
tropic patterns for snow, white sand, and stratus clouds. Stratus clouds were 
found to be more anisotropic than the other two surfaces. These clouds were sig- 
nificantly more anisotropic in the backscatter direction than in the forward scatter 
direction. Snow had higher reflectances in the forward scatter region, most nota- 
bly in the specular reflection region while white sand was more of a back-reflect- 
ing surface. For the three surface types, an increase in anisotropy accompanied an 
increase in incident angle. 

Ruff et al. [9] studied the cloud anisotropic properties from TIROS IV meteo- 
rological satellite data. This study showed clouds are primarily forward reflect- 
ing. Again, a marked increase in anisotropy was observed as the solar zenith 
angle increased. 

The deployment of a five-channel scanning radiometer aboard the NIMBUS 3 
spacecraft presented the first opportunity to study the earth's radiation budget 
with high resolution data [2]. NIMBUS 3 was launched on April 14, 1969 into a 
retrograde, sunsynchronous, nearly circular polar orbit. Equator crossing times 
were 1130 (northbound) and 2330 (southbound). In order to evaluate the satellite 
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data, a set of empirical reflectance models obtained from various sources, includ- 
ing Salomonson and Marlatt [8] and Ruff et al. [9] results, was applied to three 
scene type classifications for three solar zenith angle ranges. The models charac- 
terized the anisotropy of clear ocean, high latitude ice and snow, and cloud and 
land (cloud-land) for O’- 35*, 35* - 60*, and 60* - 80* solar zenith angle ranges. The 
cloud-land model exhibited limb-brightening in the backward scatter region(60* < 
C< 80*) and was nearly isotropic in the same region for 35’ < £ < 60*. The cloud- 
free ocean showed an increase in anisotropy with increasing incident angle. 

2-3.2. BDRF from NIMBUS 7 ERB Experiment 

Bidirectional reflectance functions using broadband observations from NIM- 
BUS 7 Earth radiation budget (ERB) data have been developed for a number of 
scene types [12-14,16]. Taylor and Stowe [12] constructed BDRFs from Nimbus 7 
data for eight uniform surface types Qand, water, snow, ice, low, middle, high 
water and ice clouds) covering a period of 61 days over regions approximately 
160 km x 160 km, called "sub- target areas". From the dataset used over the time 
period covered for their study, these authors determined that only 3% of the bins 
were unsampled and thus required interpolation. Results included anisotropic 
patterns, SWR standard deviations, and relative dispersion, which measures the 
variability of radiance within the bin. All surfaces studied exhibited an increase 
in specularity with increasing solar zenith angles, similarity in the anisotropy of 
high water and ice clouds, high backscatter for land surfaces for large solar zenith 
angles, and limb-brightening for water surfaces. Figure 2-3 [12] show the bidirec- 
tional reflectance pattern for ocean surface for SZA range 0*-26* (Fig. 2-3a) on the 
left portion of the contour plot, and for SZA 26*-37* (Fig. 2-3b) on the right half. 
These plots illustrate the presence of the sun glint region and its shift towards the 
limb as the sun moves. 

The anisotropic modeling of clouds [13] and deserts [14] utilized a reflec- 
tance model consisting of the sum and product terms of the cosines of the solar 
and viewing zenith angles, thus establishing reciprocity between these angles. 
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Figure 2-4. Staylor correlation of bidirectional reflectance with scattering angle 
for four cloud types. 
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Bidirectional reflectance patterns as * <■ . , , , 

r s a function of the scattering angle for four 

cloud types (low water, middle water , . ,. . . 

/ high water, high ice clouds) exhibit strong 

forward scatter peaks (Fig. 2-4). The correlation of the bidirectional reflectance 

w„h the scattering angle for the four cloud categories do not show significant dif- 
ferences in the shape of correlation. 

Staylor and Suttles [14] and Staylor [24] applied an analytic fit to the NIM- 
BUS 7 radiance observations to generate BDRF for desert surfaces under a dear 
sky condition. The target sites selected were the Sahara desert, the Gibson desert, 
and Saudi desert. The Sahara desert has little vegetation and moisture and has 
expansive sand dunes and sand seas. The Gibson desert, located in Western Aus- 
tralia is characterized by a mix of sand dunes, rock outcropping and arid steppe 
vegetation. The Saudi desert site is characterized by sand dunes and sand seas 
with no vegetation. Bidirectional reflectance patterns for the Gibson and Saudi 
deserts (Fig. 2-5) indicate that this surface type is primarily a back-reflector. 
Results also show that deserts exhibit varying degrees of anisotropy (e.g. Saudi 
desert is more nearly-isotropic than either Gibson or Sahara deserts). 

2-4. ERBE Operational BDRFs 

Suttles et al. [16] analyzed the NIMBUS 7 and GOES (November 1978) instru- 
ment data. The geosynchronous spacecraft GOES was positioned at 75 *W and 
made observations between 1351V and 151V longitudinally and 60 *N and 60*S 
latitudinally. The GOES radiometers detected radiances in the narrowband spec- 
tral intervals. These narrowband measurements were converted to broadband 
radiances and combined with the NIMBUS 7 ERB measurements to generate 
BDRFs for 12 homogeneous Earth and cloud surfaces. A model exists for clear 
(0-5%), partly cloud (5-50%), mostly cloudy (50-95%), and overcast (95-100%) 
scenes over ocean, land, coast, snow and desert. A single overcast composite 
model is used that incorporates the features of the overcast over ocean and over- 
cast over land scenes. The complete set of ERBE scene types are tabulated on 
Table 2-1 . The surface type is first determined by referring to a static geographical 
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Figure 2-5. Staylor and Suttles bidirectional reflectance patterns for clear desert. 



Table 2-1. ERBE scene types 



Scene 

Cloud Cover (%) 

1 

Clear over Ocean 

0-5 

2 

Clear over Land 

0-5 

3 

Clear over Snow 

0-5 

4 

Clear over Desert 

0-5 

5 

Clear over Coast 

0-5 

6 

PC over Ocean* 

5-50 

7 

PC over Land 

5-50 

8 

PC over Coast 

5-50 

9 

MC over Ocean** 

50-95 

10 

MC over Land 

50-95 

11 ” 

MC over Coast 

50-95 

12 

Overcast 

95-100 


* PC - Partly Qoudy 
•* MC - Mostly Qoudy 



map while the cloud cover category was identified by a Maximum Likelihood 
Estimate (MLE) technique [25]. The MLE technique uses a bispectral cloud frac- 
tion identification technique that uses simultaneously measured shortwave and 
longwave radiances. Figure 2-6 illustrates, for a given geographical surface type 
and viewing geometry, an example of a bispectral radiance histogram that groups 
measurement pairs into four cloud categories. Associated with each cloud class is 
a mean radiance and other a priori information and statistics that are used to eval- 
uate the probability of a measured pair belonging to one of the four cloud types. 
The cloud category which has the highest probability of occurrence is the cloud 
class selected. The bidirectional parameters (BDRFs, standard deviations of mean 
shortwave radiances, and shortwave-longwave radiance correlation coefficients) 
are tabulated in discrete angular bins, defined for angular ranges of 0 and to 
provide discrete values for modeling and presentation. Table 2-2 lists the angular 
bins used by Suttles et al. {16]. Sparsely sampled or unsampled angular bins were 
filled with models estimated by applying reciprocity, interpolation or extrapola- 
tion techniques. 

With the exception of snow scene type, the anisotropic patterns of the ERBE 
scene types show notable limb-brightening at high solar zenith angles. Ocean 
surfaces demonstrate significant specular reflection in the region of forward scat- 
tering for £ = 0 - 45* . The models increase in anisotropy with increasing viewing 
zenith angles. Clear snow is limb-darkened for £ < 53*. 

These angular models are presently being used in the inversion of satellite radi- 
ance measurements from ERBE satellites to fluxes at TOA. 

2-5. Dlhopolsky BDRF for Clear Ocean Scene 

Dlhopolsky [20] investigated the effects of higher angular resolution on the 
instantaneous albedos for clear sky over ocean scenes. The effects of ocean sur- 
face roughness on BDRF due to the effects of wind speed are also discussed in 
this study. Dlhopolsky has found that the reflected radiation is higher for a 
higher wind speed since the wave slope causes the incident angle to be larger. 



Meanwhile, a calm surface reflects a significant amount of energy in the direction 
of forward scatter. In formulating BDRFs from ERBS data, she processed ERBS 
clear sky data for the months of April 1985 through November 1985 to generate a 
reflectance matrix sorted into 5\5'x5* discrete angular bins. This resolution is a 
significant refinement compared to the first solar zenith angle bin used by Suttles 
et al.[10] which covered O' - 26' incident zenith angles. The Dlhopolsky study 
determined that the NIMBUS 7 angular bin size was not adequate for estimating 
the anisotropy in this bin range. The BDRFs for the first ERBE solar zenith angle 
bin in particular, which ranges from 0'-26‘, generally overestimated the instanta- 
neous albedos. Reflectances are also better represented in the refined angular 
bins, especially in the specular directions. Differences in BDRFs between 
Dlhopolsky and ERBE models are on the order of 10% for nonspecular directions 
and increases significantly in the specular direction for solar zenith angles less 
than 35*. 
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igure 2-6. Maximum likelihood estimate bispectral histogram for clear lan 
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Table 2-2. 

Angular bin definitions 


Bin No. 

Solar Zenith 

Viewing Zenith Relative Azimuth 

Angle 

Angle 

Angle 

1 

0.0 - 25.84 

0-15 

0-9 

2 

25.84 - 36.87 

15-27 

9-30 

3 

36.87-45.57 

27-39 

30-60 

4 

45.57 - 53.13 

39-51 

60-90 

5 

53.13 - 60.00 

51-63 

90 - 120 

6 

60.00 - 66.42 

63-75 

120-150 

7 

66.42 - 72.54 

75-90 

150-171 

8 

72.54 - 78.46 


171 - 180 

9 

78.46 - 84.26 



10 

84.26 - 90.00 





Chapter 3 


ANALYTICAL BIDIRECTIONAL REFLECTANCE FUNCTION 


The ERBE operational models tabulated by Suttles et al. [16] can be closely 
approximated by applying an analytic fit to these tabulation to obtain an empiri- 
cal form for the bidirectional reflectance. An analytic form of the bidirectional 
reflectance for clear and partly cloudy over ocean scenes is discussed, followed by 
an analytic model for the land, snow, desert, mostly cloudy over ocean and over- 
cast ERBE scene types. The forms of these expressions are based on theoretical 
considerations. 

3-1. Analytic Form of BDRF for Clear and Partly Cloudy over Ocean 


The bidirectional reflectance for clear and partly cloudy over ocean can be 
expressed in the following empirical form [26] 


(®> 0 — Cj + 


C 2 ( 1 + cos 2 y) 
(uu n ) c ’ 


C 4 (C 5 -1) 
(«m 0 ) 1 - 5 (C 5 -cosa) 2 


(3.1) 


where u = cos0,u o = cos £, y is the scattering angle, i.e. the angle through which 
the ray is turned as it is reflected, and a is the angle from the line of specular 
reflection. These angles are defined by cosy = vv 0 cos 4>-uu 0 and 

cosa = vv q cos4> + uu Q , where v = sin0 and v 0 = sin£. The second term on the 
right-hand side accounts for Rayleigh scattering from the atmosphere and for 
atmospheric absorption. The Rayleigh phase function is given by 1 + cos y and 
the parameter C 2 is associated with the Rayleigh optical depth. The parameter 
C 3 accounts for the atmospheric absorption and is affected by the presence of 
aerosols. This form for atmospheric absorption is used instead of an exponential 
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term because it is analytically tractable. The last term on the right-hand side 
accounts for the specular reflection from the ocean surface. When radiation is 
scattered from a smooth surface, such as the calm ocean, there is a sharp peak in 
the forward scatter direction. The form of the specular term was determined by 
fitting the ERBE operational models choosing an even function in terms of the 
specular angle, a, which leads to the term (C 5 - cos a) . With high winds, the 
waves over the ocean will cause the forward scattering peak to broaden while for 
calm conditions, the ocean surface will be flat thus causing the reflection to be 
sharply peaked. The first term on the right-hand side is associated with surface 
albedo and accounts for other diffuse scattering processes. The values of the C i 
coefficients are tabulated in Table 3-1. Although the coefficients tabulated on 
Table 3-1 describe the mean model since they were based on mean radiances, the 
C ( . coefficients will vary because of the variations that exist in the atmosphere and 
the underlying surface. The diffuse part of radiation varies, depending upon the 
sea state, particles in suspension, and atmospheric turbidity. The Rayleigh scat- 
tering term is affected by variability of water vapor in the air, which causes 
changes in the absorption of radiation. The sea state, which is influenced by wind 
conditions, determines the width of the dispersion of radiances about the forward 
scattering peak, thus influencing changes in the specular term [27]. Because the 
terms in the empirical formula are expressed as sum and/or product of u and 
u Q , this form satisfies the principle of reciprocity. The expression works best for 
uu Q >0.1; however, in order to normalize the expression, it will be used over the 
full ranges of u and u Q . 

The albedo can be computed from Eq. (3.1) by integrating it over the upwelling 
hemisphere (weighted by u ). The BDRF integrates to n by the normalization con- 
dition, giving 


~ -r i 3-«n 3u2- 1 

a 0 (C) = c i +C 2 u o s | 2 ~cl + 4 _ C 


3 ) 


(3.2) 
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Table 3-1 . Model coefficients for clear and partly cloudy over ocean 


Scene Type 

Ci 

Cj 

C3 

Q 


D 

l 

Clear Ocean 

0.010 

0.023 

0.800 

0.0056 

1.060 

0.011 

2 

Clear Ocean 

0.005 

0.027 

0.900 

0.0080 

1.100 

0.016 

3 

PC over Ocean 

0.040 

0.047 

0.577 

0.0080 

1.157 

0.016 


1 Applied analytic fit to ERBE operational models 

2 Tuned to alongtrack data (Dlhopolsky model) 

3 PC - Partly Cloudy 
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The last term, which is due to the forward reflection from the ocean surface, is an 
approximation to the integral, and is independent of C 4 . This is because although 
the radiance is dispersed around the forward scatter peak by an amount depend- 
ing on C 4 , the total contribution to the reflected flux is the same. The computa- 
tion of D is discussed in Appendix A. 

Another factor that induces changes in the BDRF is cloud contamination of the 
scene. ERBE scene identification algorithms classify a scene as clear if the amount 
of cloud coverage varies between 0 and 5%. In order to account for cloud contam- 
ination in clear ocean scenes, the bidirectional reflectance is taken to be [26] 


r (model) = C 6 r c + ( 1 - C 6 ) r Q (d, (/>, 0 (3.3) 

where C 6 is the effective cloud amount within the field of view. The coefficient r c 
is the reflectance for cloud, and r 0 (6,<p,Q is given in Eq. (3.1). C 6 varies from 0 to 
5% for clear ocean. The cloud reflectance used is for middle altitude water (MW) 
clouds, computed by Staylor[13]. The correlation of the bidirectional reflectance 
with the scattering angle was shown on Figure 2-4 for variable cloud heights. The 
corresponding cloud albedo has been evaluated by Green and Smith [28] and is 
given by 


a( 0 = IYqu ' 1 +2V 1 uZ~ 1 /(u 0 , N) 


(3.4) 


where Y 0 = 0.005, Yj = 1.346 and N= 1.577 for MW clouds and 


The corresponding albedo for clear ocean is computed as [26] 


(3.5) 


a (model) = C 6 r £ + (1-C 6 )a o (0 


(3.6) 
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Because some of the NIMBUS 7 shortwave data was missing or considered ques- 
tionable due to inadequate sampling, various methods to fill in these angular bins 
were implemented and flagged. Data restrictions were imposed in fitting the 
NIMBUS 7 tabulations for the development of the analytic model. The flagged 
bins were not utilized as well as bins in which uu 0 <0.1. Solar zenith angle bins in 
which the scene type is questionable were also eliminated (e.g. u 0 < 0.3 might 
either be ice or snow and not ocean surface). 

3-2. Analytic Form of BDRF for Land, Snow, Desert, Mostly Cloudy 

over Ocean and Overcast Scenes 

The bidirectional reflectance for land, snow, desert, mostly cloudy over ocean, 
and overcast scenes is given by [26] 


r (model) = d>r Ray + (mode/)] 0-7) 

where r^y is the bidirectional reflectance due to Rayleigh scattering (second 
term in Eq. (3.1) and y V, the azimuthal mean reflectance, is given by 

V = (38) 

where Ar = terbe - *ray. 1116 ERBE operational bidirectional reflectance, r ERBE/ is 
computed from the tabulation of Suttles et al.[16]. For scenes with cloud cover, 
Rayleigh scattering is reduced due to the decreased amount of atmosphere above 
the reflecting surface. For a completely overcast scene (i.e. cloud amount > 95%), 
the mean cloud tops are at 680 mb, thus 2/3 of the Rayleigh model as determined 
over the ocean is used. As the cloud amount decreases, the amount of Rayleigh 
scattering increases. The weighting factor 0 ) describes the reduction in Rayleigh 
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scattering effects due to increased cloudiness. The azimuthal mean reflectance 'F 
can be expressed in terms of the viewing zenith and solar zenith angles by [13] 

Y = A+Bx 2 (3.9) 


where Y = v f*uu 0 and x - — • Staylor [13] discusses the significance of the A 

u + u 0 

and B regression coefficients. Figures 3-1 (a,b) depict the correlation of the direc- 
tional reflectance with the zenith angles for mostly cloudy over ocean and overcast 
scenes, respectively. The regression coefficients are very close to 1 .0 for both cases 
As in Staylor's desert and cloud bidirectional reflectance models, the combination 
of the sum and products of u and u 0 satisfy reciprocity. The last term of Eq. (3.1) 
is computed from the following expression [13]: 


(model) = 


1 + KJJ3 + cosy) 


1 +A:[c 2 -2C«u 1> + (u^) 2 + i(w o ) 2 ] 


(3.10) 


where G determines where the model h^s its minimum value and K is the model 
amplitude. 

As in the clear and partly cloudy over ocean cases, restrictions on the NIM- 
BUS 7 data were imposed to exclude questionable scene types, angular bins that 
have no sampling or sparse population, and angular bins whose ADMs were cal- 
culated using any one of the interpolation techniques discussed above. 

The model coefficients for these scene types, as well as the Rayleigh weight- 
ing factors, are tabulated on Table 3-2. 

Similarly, the model albedo a (model) is computed using the relation [26] 


a (model) = coa Ray + Aa 


(3.11) 
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where the albedo contribution of Rayleigh scattering is the second term in Eq. (3.2) 
and Aa is expressed as 

r 2 

A a = 2— + 2 Bu n 1 + u - 2u In ( 1 + u Q ) + 2 ulnu - - — — 

u o o o u ° ° 1+u 

o L o 

The BDRF is finally computed by dividing the bidirectional reflectance by the 
albedo. 

The model albedos, computed from Eq. (3.11), are plotted in Fig. 3-2 as a 
function of solar zenith angle. All scenes exhibit an increase of model albedo with 
increasing solar zenith angle, although clear snow does not increase significantly 
as do the other scenes. Clear land is not a highly reflective surface due to the pres- 
ence of vegetation. Clear desert is a brighter surface since it is generally com- 
posed of light sand and little or no vegetation. These model albedos are in good 
agreement with those determined for the ERBE operational directional models. 


( 3 . 12 ) 
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Table 3-2. Model coefficients for land, snow, desert, mostly cloudy 
over ocean, and overcast scenes 


Scene Type 

A 

B 

G 

K 

CO 

Clear Land 

0.002 

0.384 

0.138 

0.650 

1.000 

Clear Snow 

0.011 

2.517 

0.675 

0.188 

1.000 

Clear Desert 

-0.003 

0.784 

0.025 

0.412 

1.000 

PC over Land* 

0.009 

0.643 

0.350 

0.875 

0.917 

MC over Ocean** 

0.025 

0.812 

0.525 

0.988 

0.758 

MC over Land 

0.030 

1.019 

0.643 

0.988 

0.758 

Overcast 

0.024 

1.530 

0.550 

0.625 

0.668 


* PC - Partly Cloudy 
** MC - Mostly Cloudy 














Chapter 4 

VALIDATION AND RESULTS 
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In order to validate the analytic BDRFs, these models are compared with the 
ERBE operational BDRFs. The effectiveness of these models in accounting for 
anisotropy at the limb when the same site is observed from a number of viewing 
directions is evaluated with alongtrack data. A brief description of this exper- 
iment is presented. Validation results for ocean under varying degrees of cloud 
cover and overcast scene are shown. The shortwave radiance standard deviations 
(SWR a) are also computed. 


4-1. Alongtrack Data 
4-1.1. Alongtrack Scan Experiment 

During limited periods in January and August 1985, the scanning radiometers 
on ERBE were rotated in azimuth so as to scan along the orbit track rather than 
crosstrack. Operation in the alongtrack mode allows for the collections of data 
applicable to radiation directionality studies. This mode, shown on Fig. 4-1, 
shows that a single site on the ground track can be viewed from a number of 
viewing zenith angles during an orbital pass. Alongtrack data are ideal for vali- 
dating the ancillary data needed to analyze radiometric measurements. However, 
the orbit geometry constraints inhibit full angular coverages. The scan mode does 
not allow for the sun within 15* of the orbit plane (i.e. no data was collected for 
relative azimuths within 15* of the principal plane) where there is strong forward 
and backward scattering. 

In order to apportion the radiance measurements, the orbit track is divided into 
16-second intervals along the ground track. An interval is approximately 108 km. 
long and corresponds to a pixel length at a viewing zenith angle of 55*. A pixel is 
assigned to the interval in which its center falls, as depicted on Fig. 4-2, regardless 
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Figure 4-1. Use of alongtrack data for multiple views of a scene. 
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Figure 4-2. Allocation of measurements to alongtrack intervals. 
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of how much of the pixel falls in adjacent intervals. Each interval has approxi- 
mately 250 pixels assigned to it. 

The scene type of the pixels are computed by the ERBE scene identification 
algorithm [1]. The interval scene type is determined by compiling the pixels 
within 10’ of the viewing zenith angle and if 80% of these pixels agree in scene, 
the scene type with the highest scene agreement is selected to be the interval 
scene. Intervals which do not meet this criterion are not used. This requirement 
assures that the interval has a uniform scene type across it. 

In a three-dimensional broken cloud field, apparent cloudiness increases with 
increasing viewing zenith angles (Fig. 4-3). With the alongtrack data, the scene 
can be identified near nadir, which is an advantage of this dataset. 

4-1.2. Validation of Fluxes with Alongtrack Data 


Alongtrack scan data is useful in determining the effectiveness of BDRFs at dif- 
ferent viewing zenith angles since a target area can be viewed from a number of 
viewing zenith angles. A target area, viewed from nadir and from off-nadir, gives 
two different radiance measurements but must be converted to the same flux by a 
correctly-modeled BDRF. 

To validate the analytic BDRF, the ERBE alongtrack radiances are converted to 
TOA flux as a function of the viewing zenith angle. For a given alongtrack inter- 
val, the fluxes are normalized to near zenith fluxes to produce flux ratios. Since 
the same area is being viewed, only at different view angles, a flux ratio equal to 
1.0 signifies correct modeling of the BDRFs. 

4-1.3. Presentation of Results 

The patterns of bidirectional reflectance functions and shortwave radiance 
standard deviations are presented on polar contour diagrams. The patterns are 
constructed for each scene type (clear, partly cloudy, mostly cloudy over ocean 
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Figure 4-3. Effect of viewing zenith angle with apparent cloudiness. 





and overcast) and solar zenith angle range defined by Suttles et al. [16]. The radial 
coordinates correspond to the satellite zenith angle 0 while the angular coordi- 
nates represent the relative azimuth angle <> between the satellite and the sun. Fig- 
ure 4-4 shows the coordinates and the angular bins of the ERBE operation models. 
The sun lies at <>=180 * such that forward scattering corresponds to 4>=0* while 
backward scattering is in the <>=180* region. Assuming that the reflection pattern 
is symmetric about the principal plane, only <>=0-180* is shown. 

Although restrictions were imposed on the ERBE tabulated models to remove 
questionable data, the BDRF and SW radiance standard deviation were deter- 
mined for all combinations of viewing and incident angles. 

The analytical BDRFs represent the mean values within the angular bins for 
n=200 realizations. These models are computed for a set of random illumination 
and viewing angles uniformly distributed over the solar zenith, viewing zenith, 
and relative azimuth angle bins. Similarly, the ERBE models that are plotted cor- 
respond to the bin mean values determined by trilinearly interpolating over the 
given random illumination and viewing angles rather than the tabulated ERBE 
BDRFs. Instrument noise, in the order of 1-2 Wm'V 1 , were taken into account in 
the radiance calculations. 

For the BDRF plots, presented in the next section, the color bar below each set 
of contour plots indicate the BDRF range (0.5 - 2.0). The left portion of a contour 
polar diagram represents the analytical BDRF while the right portion represents 
the ERBE or Dlhopolsky mean BDRF. The angular bin BDRFs for each realization 
is converted into radiance from which the radiance standard deviations are com- 
puted. Meanwhile, SWR o plots show two color bars; the left color bar gives the 
difference (Analytic - ERBE/or Dlhopolsky (clear ocean only)) in SWR a while the 
color bar on the right gives the range of the absolute SWR a. The left and right 
portions of the SWR o contour plots correspond to the respective color bars 

described above. 

Results of the validation with alongtrack data are presented as line plots of the 
flux ratios versus the viewing zenith angle for various solar zenith angle bin 
ranges. Flux ratios computed from the analytic and ERBE models are presented. 
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4-2. Clear Ocean 

4-2.1. Clear Ocean (Untuned Analytic Model) 

The model coefficients are C 2 = 0.010, C 2 = 0.023, C 3 = 0.800, C 4 = 0.0056 and 
C 5 = 1.06 (Table 3-1). In order to account for the variabilities about the mean of the 
coefficients brought upon by influences of the atmosphere and surface, statistical 
properties for the coefficients are defined. The diffuse part Cj is assumed to be 
uniformly distributed between 0.005 and 0.015. C 2 and C 4 are constant. C 3 and 
C 5 were tuned to get standard deviations that will give a reasonable match to the 
ERBE operational BDRF. C 3 has a truncated normal distribution with a mean of 
0.8 and a nominal a of 0.075. Values of C 3 that exceed 0.9 or are less than 0.7 are 
excluded. C 5 has a truncated normal distribution with a mean of 1.06 and a nomi- 
nal o of 0.020. Values of C 5 that are less than 1.01 or are greater that 1.11 are 
rejected. Figures 4-5 (a-f) depict the mean BDRF for clear ocean for six solar 
zenith angle bins for model coefficients derived from a fit to Suttles [16] tabula- 
tion. For solar zenith angle range 0-26* and $ = 0, the largest BDRF occurs at 
0 = 25 * for ERBE while the analytical BDRF is a maximum in the area of 0 = 10*. At 
small viewing zenith angles (0 < 30*), BDRF decreases as the £ increases. For large 
0, however, BDRF increases significantly with increasing solar zenith angles. All C, 
ranges exhibit an increase in anisotropy towards the limb (limb-brightening) 
except in the near zenith of the sun (specular region). This is attributed to atmo- 
spheric scatter towards the limb over a dark ocean surface. Forward scattering is 
more prominent than backward scattering. The surface is most nearly isotropic 
at $=90*. Since the BDRF values are normalized, as limb-brightening becomes 
more apparent, other viewing angles exhibit a decrease in anisotropy to compen- 
sate for the brightening. The sunglint area shifts towards the limb as solar zenith 
angle increases. The analytic and ERBE BDRFs are in good agreement. The bias 
is in the order of -0.048 and the RMS is equal to 0.118. 

Figures 4.6 (a-f) depicts the SWR o computed from the analytic model. Both 
ERBE and the analytic models do not show a trend in o from one SZA bin to 



another. However, as anisotropy increases, so does the o. The largest difference 

occurs in the forward scatter direction towards the limb for £ > 37 . For £< 37 , the 

differences in the specular region is in the order of 4-10 Wm ^sr Generally, for 

-2 *1 

the backscatter region and <j> > 35’, the difference is in the order of + 2 Wm sr . 

The effects of wind speed are depicted on Figs. 4-7 (a-b) for two solar zenith 
angle bins. The variability of wind affects how the forward scatter peak behaves. 
The amount of wind is imposed on the model by varying the value of C 5 in Eq. 
(3.1). For low wind speed, C 5 is set to 1 .06 while for high wind speed, C 5 is set to 
1.15. For very low wind, the sea surface is almost mirror-like, giving a sharp 
reflectance peak. For high winds, the sea surface is rough, thus making the for- 
ward scatter peak broader. Aside from differences in the behavior of the specular 
region, the wind does not appear to influence the BDRF pattern towards the limb 

in the backscatter direction nor near zenith. 

A comparison of the alongtrack flux ratios (Figs. 4.8 (a-h)) for ERBE, 
Dlhopolsky, and the analytic models, show that the Dlhopolsky model produces 
flux ratios that are more nearly constant towards the limb than either of the other 
two models. This is especially true in the backward scatter direction. The ERBE 
models account for anisotropy in the forward scatter direction at £ > 53 . The ana- 
lytic model flux ratios show significant growth towards the limb. For £ < 45 , the 
error at 0 = 50* is in the order of 15% and increases to 40% at 0 = 70*. The error 
increases as the £ increases. 

4-2.2. Dlhopolsky vs. ERBE Clear Ocean Model 

Figures 4.9 (a-f) illustrate the BDRF comparison between Dlhopolsky and ERBE 
for more refined solar zenith angle bins (A£ = 5"). Angular bins that were not sam- 
pled or deemed questionable by Dlhopolsky were not included in the plots, as 
shown by the black regions. The missing data includes angular bins that contain 
specular points, which were rejected by the Maximum Likelihood Estimate tech- 
nique of scene identification procedure. For SZA < 25*, the ERBE models are 
grouped into a single solar zenith angle bin. The ERBE model does not account 
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accurately for the specular peak in this range. The peak occurs at 0 = 25* for ERBE 
in 0*<£< 5* while Dlhopolsk/s model shows a shift of the specular peak with 
increasing £ as expected. For this range, the ERBE model will then overestimate 
the radiative fluxes. For £ > 25*, both models compare well although the ERBE 
BDRF in the forward scatter region has a larger magnitude. A calculation of the 
differences between the two models showed that the largest biases occurred for 
£ > 60*. The ERBE model is significantly less anisotropic in these solar angles than 
the Dhlopolsky model. 

4-2.3. Clear Ocean (Tuned) 

Using a least squares error method, the analytic model coefficients were tuned 
with alongtrack data using fluxes that were converted from observed radiances 
using the Dlhopolsky BDRF. The new coefficients are: Cj = 0.005, C2 = 0.027, C3 = 
0.900, C 4 =0.008 and C 5 = 1.10. The corresponding statistical properties are as fol- 
lows: a nominal a = 0.045 and 0.002 for C3 and C5, respectively. Note that 
because the specular term coefficient C 4 increases from 0.0056 to 0.008, the 
approximation to the integral in the albedo specular term also is changed from 
D = 0.011 to D = 0.0159. 

Figures 4-10 (a-f) compare the BDRF mean for the analytic model and Dlhopol- 
sky model for the angular bins for which Dlhopolsky had adequate sampling. 
Angular bins of 5*x5*x5* are used for presentation. Unlike the ERBE model for £ < 
25*, the region of specular reflection for the analytic model shifts accordingly with 
sun angle as Dlhopolsky model does. The Dlhopolsky model is more anisotropic 
in this region. Both exhibit limb-brightening at higher solar angles although the 
analytic model is more anisotropic in this region. At near zenith, both models are 
in good agreement. The bias and RMS are 0.005 and 0.110, respectively. 

The flux ratios (Figs. 4-ll(a-h» show that the tuned BDRF better accounts for 
anisotropy at the limb especially for C < 46’ and £> 60*. 



$t***l-- *>■ 
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4-3. Overcast 

The mean model coefficients for overcast scenes are shown on Table 3-2. The 
model coefficients are A = 0.024, B = 1.530, G= 0.550, and K - 0.625. The Rayleig 
weighting factor is 2/3 of the clear ocean Rayleigh model, assuming that the 
mean cloud tops are at 680 mb. For this study, there was no attempt to discrimi- 
nate the overcast models by cloud optical thickness, cloud height, or cloud liquid 

water content. 

The correlation of the scattering angle with the analytic BDRF is depicted in 
Fig. 4-12. Although the clouds are not classified by liquid water content as Stay- 
lor's categories [13], a comparison with his results shows that the overcast model 
falls into the low water CLW) category. The scatter diagram shows a strong for- 
ward scatter peak, minimum BDRF at 90* < y < 120‘, and a leveling off to a BDRF 
value of 1.0 at the scattering angles greater than 120’. These results agree with 

Staylor's reflectance correlations. 

The mean BDRFs are illustrated on Figs. 4-13 (a-f) for C < 66 . For £ < 37 , clouds 
are limb-darkened in the backscatter direction while for 37* < C < 46’, clouds are 
almost isotropic with 0 < 60 ‘. In higher solar zenith angles, clouds become more 
anisotropic and limb-brightening features, especially in the forward scatter 
region, become more discernible. The bias is in the order of -0.03 and the RMS is 

0.054. 

The SWR o are shown on the right portion of the contour plots on Figs. 4-14 (a- 
f). The a decrease with increasing solar zenith angles. Addifionally, as tthe 
instrument scans towards the horizon, a decreases. The left half of the contour 
plots show the differences in SWR o between the ERBE and analytic models. For 
; < 26*, the largest differences occur at the limb while for large incident angles, 
differences in the order of 25% are observed in the specular region. 

The flux ratios shown on Figs. 4-15 (a-h) show that the BDRF generally effec- 
tively accounts for anisotropy to the limb. An improvement is observed for t < 37‘ 

(Figs. 4-15 a,b). 
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4-4. Partly Cloudy over Ocean 

The model coefficients, as tabulated on Table 3-2, are Cj = 0.040, C2 = 0.0471, 

C3 = 0.577, C4 = 0.008 and C5 = 1.157. Figures 4-16 (a-f) compare the analytic and 
ERBE BDRFs. The peak of specular reflection is evident even in high sun angles. 
This specular peak shifts towards the horizon as the solar zenith angles increase. 
At this higher sun angles, weaker limb-brightening in the backscatter direction 
emerges. Limb-brightening in the forward scatter direction is very prominent 
and broad. Reflectance is nearly isotropic at the large viewing angles for $ ** 90* 
and becomes even less anisotropic near nadir. The bias is equal to 0.005 and the 
RMS = 0.110. 

The flux ratios illustrated on Figs. 4-17 (a-h) show that the tuned BDRF matches 
the ERBE and shows an improvement at £ > 60* in the forward scatter direction. 

4-5. Mostly Cloudy over Ocean 


For mostly cloudy over ocean, the mean model coefficients are A = 0.025, 

B = 0.812, G = 0.525, and K = 0.988. A comparison of the ERBE and analytic BDRF 
is shown on Figs. 4-18 (a-f). For £ < 37*, both models are nearly isotropic at 0 < 60 *. 
For £ > 37*, limb-brightening is more prominent for the analytic than the ERBE 
functions. Both models exhibit limb-brightening for £ > 53* and the forward scat- 
tering peak is more pronounced. Near nadir, less anisotropy is observed. 

The flux ratios (Figs. 4-19 (a-h)) show that the analytic BDRF is more effective 
than the ERBE models in accounting for anisotropy for £ < 53*. Flux ratios using 
the ERBE operational models do not show significant albedo growth in the back- 
scatter region as do the analytic models for £ > 53. 


ANALYTIC F.RBP. ANALYTIC „ l-RIIK ANALYTIC 



Figure 4-5. Comparison of analytic and ERBE BDRFs using nominal model 
coefficients for clear ocean (Ci = 0.01, C 2 = 0.023, C 3 = 0.800, 

C 4 = 0.0056, C 5 = 1.060). 








Figure 4-6. Comparison of shortwave radiance standard deviations from 
analytic and ERBE BDRFs for clear ocean 
(C] = 0.01, C 2 = 0.023, C 3 = 0.800,C 4 = 0.0056, C 5 = 1.060). 
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Figure 4-7. Wind speed effects on clear ocean BDRF . 
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Figure 4-8. Ratio of alongtrack SW flux using the nominal model 

coefficients, ERBE BDRF and Dlhopolsky BDRF for clear 










Figure 4-9. Comparison of ERBE and Dlhopolsky BDRF for clear ocean. 
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Figure 4-11. Ratio of alongtrack SW flux using the alongtrack-tuned model 
coefficients, ERBE BDRF and Dlhopolsky BDRF for clear ocean. 
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Figure 4-12. Correlation of analytic BDRF for overcast scene and scattering 
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Figure 4-15. Ratio of alongtrack SW flux using the analytic and ERBE UDRF 
for overcast scene. 









Figure 4-16. Comparison of analytic and ERBE BDRF for partly cloudy / 
(C, = 0.040, C 2 = 0.047, C 3 = 0.577, C 4 = 0.008, C 5 = 1.157). 
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Figure 4-17. Ratio of alongtrack SW flux for partly cloudy 











Figure 4-18. Comparison of analytic and ERBE BDRFs for MC / ocean 
(A = 0.025, B = 0.812, G = 0.5250, K = 0.988). 
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Figure 4-19. Ratio of alongtrack SW flux using the analytic and ERBE BDRF 
for mostly cloudy over ocean. 
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Chapter 5 
CONCLUSIONS 


In this study, an analytic expression for the bidirectional reflectance func- 
tion is formulated by applying an analytic fit to the ERBE operational models. The 
analytic BDRF is based on theoretical considerations. The analytic BDRF is a func- 
tion of viewing geometry and scene type. For each scene type, a single set of 
model parameters is required for application to any combination of viewing 
angles. The form of this expression satisfies the principle of reciprocity and is 
smooth in terms of the directional angles. Results are presented for four ERBE 
scene types namely clear ocean, partly cloudy over ocean, mostly cloudy over 
ocean, and overcast scenes. The analytic BDRF patterns match the ERBE opera- 
tional BDRFs. The analytic functions closely modeled the reflectances in the for- 
ward scatter direction but in the backscatter direction, the analytic models were 
slightly more limb-brightened than the ERBE operational models. 

Because these model coefficients are based upon mean radiances, the model 
coefficients have nominal values. However, due to the influences of variabilities 
in the surface and the atmosphere, the coefficients will vary. Statistical properties 
of the model coefficients were determined from which shortwave radiance statis- 
tics were computed and compared with the ERBE models. For clear ocean, the 
SWR a were of the order +- 2 Wm‘ 2 sr _1 of the ERBE SWR a. This difference 
becomes more significant in the specular region in the forward scatter direction. 
For overcast case, at solar zenith angles less than 46*, the largest difference in SWR 
a occurs at the limb for azimuthal angles between 90* and 120*. At the specular 
regions in the forward scatter direction, the difference is significant especially at 
higher solar zenith angles. 

The effectiveness of the analytic BDRF in accounting for the anisotropy at any 
viewing angle is evaluated with the ERBE alongtrack data. A given target area 
viewed from any direction must yield the same flux. The alongtrack data is ideal 
for validating the BDRFs since the experiment is designed to observe a single site 
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from multiple viewing angles. Because the ERBE clear ocean BDRF resulted in 
more significant flux increase from nadir to limb than the Dlhopolsky BDRF for 
clear ocean, the latter was used to tune the model coefficients. Although albedo 
growth is still evident for all models (analytic, ERBE, Dlhopolsky) studied, the 
tuned BDRF showed an improvement in accounting for anisotropy at the limb 
over the ERBE models. For other scenes presented, the flux ratio using the ana- 
lytic BDRF showed some improvement and /or matched the results of the ERBE 
operational models, especially for solar zenith angles less than 60 *. 

The results of this study can be used for validation and interpretation of data 
from future earth radiation budget missions such as CERES. 
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APPENDIX A. 

APPROXIMATION OF SPECULAR ALBEDO FOR CLEAR AND PARTLY 

CLOUDY OVER OCEAN 
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A-l. Specular Term A 2 = J Q — vs. cos 6 69 

(a) SZA = 0* 

(b) SZA = 45* 

(c) SZA = 60* 

A-2. Specular Term 2R$cos6 \s. cos 6 for C 4 - 0.0056 70 

(a) SZA = 0' 

(b) SZA = 45* 

(c) SZA = 60’ 

A-3. Specular Term 2 R s cos 6 vs. cos 0 for C 4 = 0.008 

(a) SZA = O’ 

(b) SZA = 45* 

(c) SZA = 60* 
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APPENDIX A. 

APPROXIMATION OF SPECULAR ALBEDO FOR CLEAR AND PARTLY 

CLOUDY OVER OCEAN 


The specular reflectance is expressed as 


r,(e,<u) 


C 4 (C S -1) 

( uu 0 ) 1>5 (C 5 - cosa) 2 


(A.l) 


where u = cos0, u 0 = cos £ , and a is the angle from the line of specular reflec- 
tion given by cos a = vv 0 cos<j> + uu 0 , where v = sin0andv o = sin£ . The direc- 
tional reflectance, Rj is obtained by eliminating the azimuthal dependence. Rj is 
defined as 



Substituting Eq. (A.l) into Eq. (A.2) gives 


(A. 2) 


R. = -A, A 


% 


12 


(A. 3) 


where 


C A C i - 1 ) 


(««„) ' 5 


and 


= \\ 


di j) 


0 (C 5 - cosa ) 2 


Since A 2 is of the form 


r ** 

J (a + b cosjc ) 2 
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(A. 4) 


where a = C$-uuq and b = -Wq, 


a U/L 

2 2 2 1-5 (A. 5) 

(a-b 2 ) 

Figures A-l (a-c) depict A2 as a function of u for three solar zenith angles, 0 = 0*, 
30*, and 60*. These curves are the same for any value of C4 since A2 is independent 
of C4. The peak in each solar zenith case occurs where u = uq . The smaller the 
value of C5, the greater the area under the curve. 

The specular model albedo is computed by numerically integrating Eq. (A.2) 
over the viewing zenith angle as 

A(d Q ) - 2j l Q R s udu (A. 6) 

Figures A-2 (a-c) and figs. A-3 (a-c) for C 4 = 0.0056 and C4 = 0.008, respectively, 
show 2R s u as a function of u for varying values of C5 for the three incident angles. 
As C4 increases, the specular term increases. The peak for each case occurs at u = 
u Q . The C 5 is associated with surface roughness in the specular term. An increase 
in this parameter can be interpreted as an increase in wind speed and a decrease 
corresponds to low or calm conditions. From this figures, a low C5 corresponds to 
high specular peak while a large C5 corresponds to a broadening of the specular 
term. 

For each solar zenith angle, albedo is computed from Eq. (A.6). D is computed 
as 


(A. 7) 


D = jjZA(e Q )u 2 Q 


where N = 3 (i.e. 0 O = 0*, 30 >0*). 


!5 


i uu 
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APPENDIX B. 

COMPUTATION OF MODEL COEFFICIENTS FOR CLEAR AND PARTLY 

CLOUDY OVER OCEAN 


Fi gure 


Page 


B-l. Specular Term R s uuq vs. cos a 
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APPENDIX B. 

COMPUTATION OF MODEL COEFFICIENTS FOR CLEAR AND PARTLY 

CLOUDY OVER OCEAN 


Proposed Model: 

C, ( 1 + cos 2 y) C 4 (C 5 -1) 

c i + ( U u 0 ) Ci + (uu 0 ) 15 (C 5 -cosa) 2 


Data Source: ERBE operational models [14] 


(B.l) 


1 . 


Compute bidirectional reflectances from tabulated ERBE BDRF R e (6,<f),Q and 
mean albedo a(Q 

r e = R e (e,4>,0a(0 (B.2) 


2. Data restrictions: 

Exclude the following data - 

(a) uuq < 0.1 (non-plane-parallel) 

(b) uq < 03 (questionable scene type) 

(c) flagged data (insufficient sampling) 

3. Reduce the specular term to a small value by editing out cos a > 0.6 

giving the diffuse surface and atmospheric reflection terms only as 


(0, <>, 0 » Cj + 


C 2 ( 1 + cos 2 y) 

(KU 0 ) C ’ 


(B.3) 
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4. Iterate on Cj and perform linear regression of 

f R-C ' 1 

l n i Li vs. In (uu 0 ) 

1 + cos y 

where R = r# until the highest regression coefficient is obtained. Check for lin- 
ear relationship. 

5. Solve for initial values of Cj , C2, and C3. where 

C2 = exp(intercept) 

C3 = -slope 

6. Check for the validity of the specular model. The specular reflectance is 

r s- r e~ r d 

Figure B-l shows the specular reflectance term as a function of cos a. 

7. Plot and regress lrs(uuO) 15 r° 5 vs. cos a for cos a > 0.6 only. Solve for 

_i , - -intercept 

(2 — — - End Cj — : 

4 slope ( slope + intercept ) slope 


8 . Recalculate rA6,<p,0 = ( r e ~ O where 


r f (d,< p,q) = 


C (C - 1) 

4 —1 and recompute Cj, C 2 , and C3. by 


(uuo) 1 - 5 ^ -cosa) 


repeating steps 4 and 5 and substituting R = rj. 



cos a 


Figure B-l. Specular Term R s uu 0 vs. cos a 
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APPENDIX C. 

COMPUTATION OF MODEL COEFFICIENTS FOR LAND, SNOW, DESERT, 
MOSTLY CLOUDY OVER OCEAN AND OVERCAST SCENES 
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APPENDIX C. 


COMPUTATION OF MODEL COEFFICIENTS FOR LAND, SNOW, DESERT, 
MC / OCEAN AND OVERCAST SCENES 


Proposed Model: 


r (model) = (or Ray + (mode/)] 


(C.l) 


where 


(model) = 


1 + K(G + cosy) 


^ 2 1 
1 +k[g 2 -2Guu 0 + (uu 0 ) +2^ vv o^ J 

Data Source: ERBE operational models [10] 


(C.2) 


1 . Compute bidirectional reflectances from tabulated ERBE BDRF R e (6,<t>,0 and 

mean albedo a(Q 

r e - R e (ft h 0 a (0 

2. Data restrictions: 

Exclude the following data - 

(a) uuq < 0.1 (non-plane-parallel) 

(b) SZA restrictions: 

clear land -uq<03 and uq > 03 
clear snow - uq > 0.7 
clear desert - uq < 0.4 

(c) flagged data (insufficient sampling) 
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3. Compute the azimuthally dependent reflectance by 

•F = ±f*Ard<|> (C.4) 

TCjU 

where 

= r t~ r ray (C5) 

r ray is the clear ocean Rayleigh reflectance term. Compute (data) 

4. Regress ^ (data) vs. ^ (model) varying G and K until maximum 

regression coefficient is obtained. 

5. To calculate A and B, 

Let AY = and x = uu 0 / (u+Uq) 

Perform linear regression of AY and x 2 . Solve for A and B from 
AY = A + Bx 2 


